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Abstract 

We derive the kinetic equations for the hopping transport that take into account electron spin 
and the possibility of double occupation. In the Ohmic regime the equations are reduced to the 
generalized Miller-Abrahams resistor network. We apply these equations to the problem of the 
magnetic moment relaxation due to the interaction with the random hyperfine fields. It is shown 
that in a wide range of parameters the relaxation rate is governed by the hops with the similar 
rates as spin precession frequency. It is demonstrated that at the large time scale spin relaxation is 
non-exponential. We argue that the non-exponential relaxation of the magnetic moment is related 
to the spin of electrons in the slow-relaxing traps. Interestingly the traps can significantly influence 
the spin relaxation in the inhnite conducting cluster at large times. 

PACS numbers: 
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I. INTRODUCTION 


In recent years investigations of the spin phenomena in hopping transport have gone 
through a vigorous revival. In particular, the increasing interest to this problem is related 
to the discovery of a strong spin-valve effect in organic semiconductor devices [l, 2|. There 
is a consensus that the conductivity in these devices is determined by hopping polaronsj^. 
Such devices show a number of spin-related phenomena, including the spin-valve effect itself 
that are not understood. These experimental results were followed by a number of theoretical 


investigations of the hopping conduction 


4 l5| including fluctuations of conductivity jb, ^ 


and most important the spin phenomena like magnetoresistance 8l4l2l| and the relaxation 


of magnetic moment 


13 


19|. 


Most of these theoretical studies are based on semi-qualitative concepts with the lack of 
solid theoretical proof. Although in some cases (for example in the case of strong Coulomb 
interaction) the self consistent theory of hopping transport does not exist, there is a way 
to make theory of hopping transport self consistent. This theory is described in Ref. 20 |. 


Up to now t 
density 


21 


re theory includes the electron spin only in the limit of the small electron 


23|. 


Here we develop the generalization of this theory to include spin and arbitrary probability 
of the site occupation (with possible double occupation). We argue that the self consistent 


io| 


theory of hopping transport is necessary in order to test the semi-qualitative concepts 
and to understand the electron and the spin transport in organic semiconductors. 

Our starting point is the general Hamiltonian that describe a system of localized sites 
with the possibility of hopping due to electron-phonon interaction. Then we derive general 
kinetic equations that describe both the charge and the spin transport in hopping media. 
We argue that these equations are a useful tool to study hopping transport phenomena. 

In the present paper we apply these kinetic equations to the problem of the spin relax¬ 
ation due to random on-site magnetic helds in a system with the positional disorder. The 
understanding of the spin relaxation is clearly important for the theory of the spin-related 
transport phenomena like spin-valve magnetoresistance. Up to now the most studied mech¬ 
anism of the spin relaxation in hopping transport is the spin-orbit interaction. The theory of 


this relaxation mechanism is discussed in Refs. 


13 


17j . In Ref.[l3| the basic understanding 


of spin relaxation for the hopping transport due to spin-orbit interaction was formulated. 
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In Ref.[l^ the results of were signihcantly expanded. It was pointed out jl^ 
exponentially broad distribution of hopping rates plays determinant role in the relaxation. 
However this distribution was considered semi-quantitatively. In Ref. it was proposed 
that in the case of hopping conduction the spin relaxation due to the spin-orbit interaction 
may be signihcantly enhanced by the exchange interaction. 

While the spin-orbit interaction is supposed to dominate the spin relaxation in non- 
organic semiconductors with hopping conductivity the situation in the organic materials may 
3e different. In organic semiconductors the spin-orbit interaction is substantially suppressed 
16l. llTI. |24| and another mechanisms like hyperhne interaction can govern the spin relaxation. 
The theoretical description of this mechanisms is different because the spin-orbit interaction 
does not affect the spin of a localized electron and manifest itself only in the spin rotation 
during the hop. On the other hand the hyperhne interaction leads to the appearance of 
random effective on-site magnetic helds that rotate electron spins even without hops. Note 
that in the organic spin-valve devices there is the another source of random helds unique 
to these case. The hnite roughness of the contacts leads to a leakage of the magnetic held 
from the ferromagnetic contacts to the organic layer 25(]. This fringe magnetic held can be 
a source of additional mechanism of the spin relaxation. In terms of theoretical description 


it is added to the hyperhne held and a 


Recently Harmon and Flatte 


18 


so rotate spins on the localized sites. 


19| proposed an interesting approach to the spin relax¬ 


ation based on the waiting time distribution. They considered both the spin-orbit and the 
hyperhne relaxation mechanisms. However their approach does not take into account the 
Pauli pri ncip e and is applicable only to the limit of the small electron concentrations. Also 
Refs. [l8|, ll9| consider only the energy disorder and the positional disorder was neglected. 

The derived kinetic equations allow us to describe spin relaxation with any site occupation 
probability. In the present work we apply them to the problem with the positional disorder 
and neglect the energy disorder. We show that even in this case there are several phenomena 
in the spin relaxation that were not discussed previously. Most important of them is perhaps 
the non-homogeneous character of the spin relaxation. In some cases the relaxation of spin 
of the electrons important for conduction is substantially different from the relaxation of the 
average spin of the system. 

The paper is organized as follows. In section [TTl we discuss the usual way of describing 
hopping conduction and its justihcation in terms of kinetic equations. In section IHII we 
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generalize approach 20(] to include electron spin and derive general form of kinetic equation. 
In section IIII Bl we linearize the kinetic equations and derive generalized version of Miller- 
Abrahams resistor network. Finally in section HVl we use the kinetic equations to describe 
the spin relaxation in the hopping conduction regime. 


II. KINETIC EQUATIONS FOR HOPPING TRANSPORT 


The usual approach to the theoretical description of the hopping transport starts with 
the introduction of the hopping rates between pairs of sites 2^. The hopping rate from site 
i to site i is dehned as 

r,,cx|4f(1) 


Here Itj oc exp(—r^j/a) is the overlap integral between sites i and j, rij is the distance 
between these sites, a is the localization radius. and are the occupation probabilities 
of sites i and j. AEij is the difference of energies of states i and j. Pph{AEij) is the part 
of hopping rate related to the number of phonons, involved in the hopping. It is equal to 
Nph + 1 when AE^j > 0, i.e., when the hopping occurs with the phonon emission and to Nph 
otherwise. Here Nph is the number of phonons with energy \AEij\. 

Then the current between sites i and j is introduced as 


Jij — ^(Xij Tj-j). 


( 2 ) 


Jij = 0 when the external electric held is absent. In a weak electric held the current follows 
the Ohm law = Uij/Rij with ehective resistor voltage Uij = Acpij — A/i^/e corresponding 
to the shift of electrochemical potential between sites i and j. The resistance Rij is dehned 
as 


Rij 


knT 


e^r 


( 0 )’ 


(3) 


where T® is the hopping rate between sites i and j without the external held, fcs is the 
Boltzmann constant. 

As a result the real system with the hopping conductivity is replaced by a network of 
classical resistors. This network can be treated with percolative methods. In this case 
one hnds the threshold resistance that allows the percolation to the macroscopic distances. 
In the case when the distribution of resistances Rij is exponentially broad, this threshold 
resistance governs the conductivity of the whole system. 
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Without Coulomb interaction this approach can be consistently derived in terms of kinetic 


equations 


20j. It is the starting point of many theoretical studies and works quite well 


in many cases. However it has some problem when the hopping conductance becomes 
dependent on the electron spin because spin is not included in the equations ([H-E]). Let us 
discuss how this scheme can be modihed in order to include electron spin and how it was 
altered in previous studies. 

The essential part of the discussed treatment is the percolation theory that is a conven¬ 
tional way to study analytically dc current in a hopping system with a broad distribution 
of hopping rates. As long as the hopping rates depend on spin and electron spins are not 
totally polarized the theory should also contain some sort of averaging over possible spin 
direction. There are at least two possible orders of this procedures: one can hrst do the 
spin averaging and then calculate the percolation parameters or calculate the percolation 
threshold hrst and then average over spin (or make some more sophisticated approach to 
relating these two procedures). Different existing studies apply different order of percolation 
and spin averaging. 

In Ref. the percolation was calculated over resistors network that considered both 
probabilities of having spin up and spin down electron on each site. These probabilities 
depend on the magnetic held when Zeeman energy is of the order of kT leading to the 


ation with rates averaged over spin directions was also im- 


magnetoresistance. The perco’ 
plicitly considered in j^| and 


29| where it was noted that electron spins should be freezed 
in order to have negative interference magnetoresistance (otherwise the averaging over spin 
directions kills the ehect). 

The another relation between percolation and spin averaging was considered in Refs, js- 


Q 


In these articles it was allowed double-occupation of the sites but only for electrons 
with antiparallel spins (more exactly: in singlet spin state). This assumption agrees with 
but the following a ppr oach is diherent. 


The approach in 


lfll | starts with momentary site occupations and momentary spin 


projections. Then the percolation is considered in terms of effective concentration of sites 
allowed for the hop of a given electron. Double occupied sites are always excluded from the 
percolation (disregarding the fact that they may be empty in following moments of time). 
Single occupied sites are included in the percolation when the spin on the site is antiparallel 
to the spin of the hopping electron. If electron spins are parallel the site is included in 
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percolation with some probability p{H) reflecting that the spin flip can occur faster than 
the hopping of an electron to a distant site. 

With p = 0 the percolation is calculated before spin averaging. Finite p{H) leads to a 
more sophisticated relation between percolation and spin averaging. The dependence of p 
on magnetic fields leads to the magnetoresistance. 

The approach different from previous two was proposed by Osaka |30|. The percolation 
in the model of the resistor network was considered. The spin-flip process was considered 
as an additional resistor connected in parallel to a normal one. 

We want to underline that the discussed approaches are not equivalent and lead to dif¬ 
ferent physical results. For example in the approach of 27| the concentration of sites that 
participate in percolation is the full concentration of sites as long as there is some probability 
for any site to be allowed for the hop. In js-lO] at least double occupied sites are always 
excluded from the percolation. 

To understand what is the correct procedure of spin introduction into the percolation 
theory let us consider the procedure of the consistent derivation of the approach ([T]). This 
procedure is described by Bottger and Bryksin j^, however we remake it here in a slightly 


different way to make its generalization easier. Note that the discussion in Ref. 


20j is 


focused on the case of low occupation numbers and it is stated that that the result for any 
occupation number is similar if the Hartree decoupling is used. We make our theory with 
Hartree expansion from the beginning.. 

We start with the hopping Hamiltonian after polaron transformation, which is well known 
in the polaron transport theory 20(]. 


H — Hq + Hhop, Hq — — Ep)al Oj -|- fn/jq + 2 j ’ 

i q ^ ^ 

Hjiop ^ tjjQjj^ 

i.,, = exp - ft-c.]} 

q 

Here a, is the electron annihilation operator on site i, bq is the annihilation operator for 
a phonon with the wave-vector q and with frequency Uq. We take into account an energy 
disorder therefore each site has a random energy tij are the overlap integrals between 
the sites, t^j are much smaller than the differences in random energies e, — £j. Mq is the 
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electron-phonon interaction matrix element Mq = Mlq. r* is the position of the site i, 
Ep = iMqp/hcjq is the polaron binding energy. 


The conventional derivation of the kinetic equation Ml] starts with the full description 
of the system (with Liouville equations if the system is classical or with Hamiltonian if the 
system is quantum). Then the system is divided into the parts that interact weakly or rarely. 
We divide the system into the set of noninteracting sites, described by the Hamiltonian Ho, 
and weak interaction of different sites, described by the Hamiltonian H^op- In order to 
describe the system with the set of kinetic equation we have to assume that the eigenstates 
of the Hamiltonian Hq are well dehned and therefore H^op is much smaller than Hq. 

In zero order over H^op fhe full density matrix of the system is diagonal in terms of 


electron hlling numbers rij. The theory of hopping conduction near Fermi level 


that we are 


interested in) is based on the Hartree-Fock decoupling for the density matrix 20||. With this 


decoupling the full density matrix can be expanded as a product of the single site density 


matrices. 


_ „[{!) . n'(2) . . 

Pni,n2,...,njv — Hrii Hn2 


(5) 


Here we used the basis of the Filing numbers. The set of the system states that has dehnite 
hlling numbers n, is the complete set of states of the system. As usual the density matrix 
has two indexes each of them corresponding to one of the states of the complete set. The 
upper indexes in round brackets correspond to a number of site. 

The single site density matrix pn\ has only two indexes n, = 0,1 that are the possible 
hlling numbers of site i. Without Hhop this matrix is dehned as: 

/« 0 
0 1 - /(*) 

where is the probability for site i to have an electron. 

n'- (i) 

Here we will use the interaction representation, therefore pnl does not depend on time 
without perturbation. In the hrst perturbation order we obtain 


p« = 


l(i) 0(i) 

Pi Pi 

l(i) 0(i) 

do Po 


( 6 ) 




1 


dt ih 


{Hhop)ij,p 




p{ij) ^ p{i)p{3) + 


(7) 


Here {Hhop)ij is the part of Hhop corresponding to the hops between sites i and j in the 
interaction representation. jg the two-site density matrix. We divide p*^uO j^to the 

product of the one-site density matrices and a small oh-diagonal (in terms of electron hlling 
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numbers) correlated part oc Hhop- The contribution of the product corresponds 

to the mean-held correction to the energy Si and is usually neglected. The time evolution of 
p(*) ig governed by the off-diagonal part 

To obtain the equation for one should write the equation for the two site density 
matrix 


dt 


ih - 


(^Hfiop)ijy ^ ^ Tr^ (^Hhop)kjy P 


(ijk) 


( 8 ) 


Note that enters equation ([7]) with the coefficient oc Hhop- The right hand side of Eq.(|8]) 
also contains Hhop- Here we neglect all powers of Hhop higher than 2. Therefore we replace 
all many-site density matrices by the products of the one-site density matrices in the right 
hand side of Eq.(IH]). As a result we obtain the expression for d*'*T(t); 






dt'. 


(9) 


To get the kinetic equation we assume that the correlated part decays much faster in 

comparison with the time scale on which the one-particle density matrix changes. Therefore 
the one-site density matrices are out of the integral in Eq.([9]). This assumption allows to 
get the hnal equation for the one-site density matrix: 


dt 




i^HflQp )ij{t)y 


{Hhoph{t')dt',p^^\t)p^^\t) 


( 10 ) 


ph 


Hhop contains not only electron but also phonon operators. Therefore we should average the 
right-hand side of this equation over phonons. Substituting Eqs.(jl]) and ([6]) into flTOl) we 
obtain the following result: 

d/w 


dt 


_ /W) _ - /(T), 


( 11 ) 


where the hopping rate in the limit Mg/hWg 1 has the form: 




|M,|" je"’'- - [(iVpt + l)i(£.: - Ej + Ihj,) + N^^Sie, -Sj- fio.,)]. 

q 5 

_ ( 12 ) 


The hoping rates in the limit of Mg / hujg > 1 are derived in 201 


It is the kinetic equation for the spinless electrons. In order to derive the dc current one 
should consider the steady state solutions of these equations dfi/dt = 0 and then introduce 

















a small electric field. As a result the equation yields the standard formula for the current 
© where h,,- = IV,,■/«(! - 

Let us discuss the applicability conditions of the kinetic equation. First of all it re¬ 
lies on the smallness of the intersite correlations In the discussed problem it is 

closely related to the Hartree decoupling. It is valid when the system is close to equi¬ 
librium 32|]. However far from equilibrium the site occupation may be correlated, therefore 


{ajttiO^ttj) 7^ {af tti) {Uj ttj). In this case the system cannot be described with a closed equa¬ 
tion for /,. Real systems, however can have other sources of these correlations (even in the 
equilibrium). The most known of them is the Coulomb interaction. The problem of the 
Coulomb interaction in hopping conductivity is rather long-standing. If the Coulomb inter¬ 
action between neighboring sites is comparable with the random energies, the correlations 
can appear even without Hhop- In that case the applicability of the kinetic equation becomes 
questionable even in the spinless case. The standard answer for this problem is that the 
kinetic equation is still applicable but the density of single electron states should be consid¬ 
ered taking into account the Coulomb interaction. It leads to the formation of the Coulomb 
gap near Fermi level. Equation (1TT|1 with the Coulomb gap was successfully used to obtain 
the well-known Efros-Shklovskii temperature dependance of conductivity 2^, however the 
other consequences of Coulomb interaction are still under discussion. 

Another important simplihcation is neglecting of all the high-order terms of Hhop/^i- It 


is known that these higher order terms may leac 


to the phenomena that are important for 
33j it was argued that intersite exchange 


the hopping conductivity. Recently in Refs. jlS, 
effects may be important for the spin relaxation and the spin transport. Another example 
is the sub-barrier scattering that is responsible for the linear negative magnetoresistance in 


semiconductors with the variable-range hopping conductivity 


28| . Although these phenom¬ 


ena are important, we believe that it is more important to understand low-order physics 
before consider these high-order phenomena. 


III. KINETIC EQUATION WITH SPIN AND POSSIBILITY OF DOUBLE OC¬ 
CUPATION 


Up to this moment our results repeat the traditional scheme 2q , at least when dc current 
is under consideration. The advantage of our approach is that it may be easily generalized 
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for the case when electrons have spins and can double occupy a site. 

Equation ffTOj) does not depend on the exact structure of the one site density matrix 
corresponding to spinless electrons. It relies only on the kinetic equation assumptions that 
are valid in general case. Therefore in order to generalize our theory we have to include spin 
into the Hamiltonian, dehne the structure of one-site density matrix in general case, and 
perform the calculations that are analogous to the derivation of Eq. lfTT]) . 

Here we consider the case when the electron spin conserves during the hopping (the 
under-barrier spin rotation is discussed in section UlI Ap . However we include small on-site 
spin Hamiltonian Hs that describes rotation of the spin over the effective local magnetic 
held Hs = ^ibg Yji HiS j, where Hj is the local effective magnetic held and % is the operator 
of spin on site i. 

H = Hq + Hs + Hhop- (13) 


Hq ^ ^ T ) Hs ^ ^ Hs^i 

i,s i i 

Here s is the spin index that can have two values “-I-” and ”. U is the Hubbard energy. The 
introduction of the term corresponds to the following model. We allow 

double occupation of the site but consider that other excited orbital states at the same site 
have very large energies. So far both electrons on a double-occupied site have the same 
coordinate wave-function and their spins should form a singlet state. Double occupation of 
a site with two electrons in the triplet spin state is not allowed. Hs^i = acts only 

on electrons on site i and conserve the hlling number of this site, i.e., it can only rotate the 
spin on a single-occupied site. 


Hhop — (14) 

ijs 

where ^ij is dehned after Eq.(j4]). The hopping part of the Hamiltonian (IT4)l conserves the 
spin. 

i'_ ,i' 

The one-site density matrix in the representation of hlling numbers now contains 

four indexes and have 16 matrix elements. We however will use another representation in 
this section. The complete set of states for one site with possibility of double occupation 
has four states that can be selected as 


|0), |+), I-), |2). 


(15) 
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Here |0) is the empty site, |2) is the double occupied site, |+) and |—) are the states of the 
site when it has one electron with spin up and down respectively. So one can write one-site 
density matrix with only two indexes, each of them can have any of four values described in 

m- 

In zeroth order over H^op we have only matrix elements of the one-site density matrix that 
conserve the hlling number. For the spinless electrons there were only two such elements. 
Now we have six of them. 

Po, pX, P+, pt, pi, pI- (16) 


Only these 6 matrix elements appear in the kinetic equation. 

Along with the density matrix elements ffT6D it is sometimes useful to consider another set 
of 6 numbers, that are linear combinations of the matrix elements (1T6|1 . in order to describe 
the state of the site. These numbers allow us to track directly the occupation number and 
the mean value of the magnetic moment of the site. The transition to these numbers was 
proposed by Bryksin (without double occupation probability) [3^ . 


/o = Po, fi = pX + p_, 

M^=pX-pZ, M^=pt + p+, (17) 

My = -ip-^ + ipt, f2 = pl 

Here /o, /i and /2 are the probabilities for a site to have 0, 1 and 2 electrons respectively. 
Ma is the mean magnetic moment of the site in the direction a. 

The generalized equation ffTOj) in these notations has the following form: 

^ w,,Si)pft>Z- ( 18 ) 

j 

Here x, y and z have 6 possible values, p^ correspond to some hlling number probabilities 
or mean projections of magnetic moment dehned in (1T7|1 . The term Sxy{i)py'^ corresponds 
to the action of the spin part of the Hamiltonian Hg- 


S 


PbQ 




(ii) 


Here eay^y is Levi-Civita symbol, pb is the Bhor magneton and g is the g-factor. It correspond 
to the precession of the local magnetic moment dWi{i)/dt = {pbg/K) [M(^))H(i)]. We also 
introduced here a renormalized local magnetic held measured in units of frequency = 
(PbP/h)HW 
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The important assumption related to Eq. ffTSD is that the on-site Hamiltonian can be 
treated independently from Hhop- It is valid when Hs5tMark/f^ -C 1, where Stuark is the 
characteristic decay time of the correlation dij. (In the opposite limit the kinetic equation 
becomes non-Markovian). We also did not include the Zeeman energy when the averaging 
over phonon states Eq. (ITUD is discussed. This assumptions are justified when Hs is small 
compared to Si — Sj and kT 

For larger magnetic helds especially when Zeeman energy becomes larger than temper¬ 
ature the kinetic equations in the present form are valid only when the site magnetization 
and the magnetic held are oriented along one axis (in that case Hs does not lead to the 
magnetization precession). In that case the site energies in the equation should include the 
Zeeman energy. If the on-site magnetic held is large and is oriented along diherent axes the 
phonon averaging can lead to more complex equations. This case is however out of the scope 
of the present work. 

The term Wxyz{'ij)py ^is the “collision integral” 


1 

W 


Tr, 


i^HfiQp^ 




(19) 


ph 


represented in the notations (ITT)) . Each of the indexes x, y and z can have 6 diherent values, 
so there are 6^ matrix elements Wxyz and their calculation is rather cumbersome. Using the 
trick described in the Appendix we derive the following set of kinetic equations: 




h) 


dt 




(i) fU) , 
0 


IU,“ 


-u 




fPfP - 


( 20 ) 


- 2w;pfPfP. 

^ = E +2Wiy/<■>/“>- 


( 21 ) 


T-U 


W-^ + W- 

V ^ 


-u 


2 

dMii'> 






dt 


+ = E ’ + M-'yMf’/f+ 


( 22 ) 




W, 


-u 

ij 


+ lU" 


-U 




df, 


(*) 


dt 




(i) fU) ^ 


IU“ 


-u 




(23) 
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Here we have introduced the hopping rates and One can note that Wij 

defined in (IT^ depends on the energy difference £* — Sj. However if one of the initial and 
the final state of the hop corresponds to the upper Hubbard band, the actual energy Si + Uh 
or Ej + Uh should be used. Therefore is the hopping rate Wij where Si is substituted 

by Si + Uh and W^^ is Wij with Sj substituted with Sj + Uh- 


A. Spin-orbit couplings 


The kinetic equations Eqs. fl20tl23p were derived with the approximation that the electron 
spin is conserved during the hop. This approximation is not sufficient when the spin-orbit 
interaction is essential for the kinetics. Although we do not want to discuss the role of 
the spin-orbit interaction in details we briefly outline the procedure of the inclusion of the 
spin-orbit interaction in the kinetic equations in this section. 

The possibility of inclusion of the spin-orbit interaction into kinetic equation for hopping 
conductivity was discussed in 2ll-l23|. where corresponding kinetic equations were derived 


in the limit of low occupation numbers. The spin-orbit interaction leads to a precession of 
the electron spin during the under-barrier motion. It is important that for a given pair of 
sites i and j the angle of precision is not random (it is the same for all hops between these 
two sites). Therefore the spin orbit interaction can be described with rotation matrixes 
If the electron on site i has momentum expectations Ma (where index a stands for 
the cartesian coordinates), then after the tunneling to site j the expectation value of the 
magnetic moment is This rotation matrix should be defined for every pair of 

sites 

To include rotation matrixes into kinetic equations (1201123^ . one should make a substitu¬ 
tion 

AifP ^ UfW (24) 

in all kinetic equations Eqs. (120112^ . On the other hand projections Ma'^ should be un¬ 
changed. 
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B. Generalized resistance network 


When the linear response of a hopping system to a small applied dc voltage is considered 
the kinetic equations can be reduced to a resistor network. In this section we show how 
this network is generalized when the electron spin and the Hubbard energy are taken into 
account. 

The reduction to the resistor network depends on the possibility to introduce quantities 
that change slowly in space. Note that probabilities ft and magnetic moments Mi differ 
signihcantly from site to site even in the equilibrium in the case of the broad distribution of 
random energies. 


In the case when the directions of magnetizations on all the sites are the same 


35j one 


can introduce chemical potentials for spin up and spin down electrons /i+ and /i_. The 
occupation probabilities should be expressed as functions of these chemical potentials. 




ksT 


P-M = 


-Ci+M. 


(^) 


_ 

J 2 — 


exp 


-2ei-Uh+IMf 

kgT 


where Zi is the statistical sum on site i 

ry ^ — Ci + pV —Ci + —2ej — [//i +/i^ + 

Zj = 1 + exp —-——-h exp —-——-h exp- 


(25) 


(26) 


hj^T ksT k]^T 

Here we assume that the magnetization of all sites is directed along z axis, therefore = 

p+iP ^ g_ Jii) 

are expressed in terms of and according to Eq. flTT)) . 

It is easy to check by the direct substitution that when chemical potentials are the same 
in all sites p_|_ = /i_ = /i the system is in the equilibrium for any on-site random energies 
and all time derivatives in the kinetic equations are zero. 

It is useful to have expressions for the spin up and the spin down currents. The spin up 
current between sites i and j is given by the formula: 


4 = -'(rj - rt), 


Here Th is the hopping rate for the spin up electrons from site i to the site j. 


r; = /t’ + 


c{i) Aj) 


(i) Aj) 


r-U AP f(j) 


(27) 


(28) 


To get this expression from the kinetic equations one should consider the time derivative 
of the probability to have an electron with the spin up on site i. Note that a double 
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occupied site contains electron with spin up and thus + f!f^. Direct calculation 

yields: 

P9) 

j 

with Eq. fl27D for the currents The expression for spin down current can be obtained 
from fl271) and fl28|) by replacing index + with —. 

The current appears when the system is placed in the electric held or when there is a 
shift of the chemical potential. Expanding the kinetic equation over the small electrostatic 
potential Acpi and the small shift of the chemical potentials Afif we obtain: 


J* = - A^^ye) 


(30) 


where A/rd = fi'l — , Acfij = ipi — ipj and is the effective resistance of the pair of sites 

ij with respect to the current of the electrons with spin up. It contains four contributions. 

r 1-1 


-Ri/,+ 


kuT 


(31) 


where Td with additional indexes A and B are the contributions to Td 





The additional upper index (0) in fl3ll) means that values T are calculated for Aipij = 0 and 

Aflij 0 . 

The physical meaning of different contributions to T becomes apparent when one considers 
large Hubbard energy U 3> ksT. Note that in order to contribute to hopping conductivity 
a site should have the energy level close to the chemical potential. In this case the sites 
that take part in conductivity are divided into two groups j^: A-sites that have e* ~ /i 
and B-sites with £i-\-U ~ /t. A-sites are practically never double-occupied so for these sites 
/2 ~ 0. For B-sites /o ~ 0, these sites never have zero electrons. For this model only one 
contribution to is important for each pair of sites. For example if the site i is of type 
A and the site j is of type B, the resistance for this pair is 


R 


V, 


ksT 


e^T 


+,(o) • 

ij,AB 
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This result agrees with 27|] when the temperature is much less than the Hubbard energy. 
We generalize this result to the case of arbitrary relation between kT and U and make it 
explicitly applicable to the case U = 0 discussed in js-lO]. We show that if the temperature 
is larger or comparable with the Hubbard energy each site plays both roles: of A-type and 
B-type. The corresponding resistances are connected in parallel. 

When magnetization is not restricted to one axis and the magnetization of sites in differ¬ 
ent parts of the sample is aligned along different axis, the introduction of slowly changing 
chemical potentials is possible only when this axis slowly changes in space. In this case the 
chemical potentials fif and /i“ should be related to the mean spin projections on the local 
axis. In this case one can use the same expressions for resistors, but should keep in mind 
that the system kinetics cannot be reduced only to the resistor network. The situation when 
different parts of the sample have different magnetization axis is unstable and leads to the 
spin relaxation. 

Finally let us note that although magnetization Mj depends (even in the equilibrium) on 
random energies e*, the relative magnetization Aii = Mi/fi is the function of the chemical 
potentials only. 


Aii = tanh 


N - 
ksT 


(32) 


C. Magnetoresistance 


A new mec 
semiconductors 


lanism of singlet magnetoresistance was recently proposed in organic 


lOl. This mechanism is based on the fact that two electrons with the 
.pin cannot occnpy the .a.e cite even .n the caee of entaU Hnhbatd enetgy. In Ref.^ 
[10] it was stated that when the spin relaxation time is longer than the hopping time, in 
order to hop the electron should hnd the site which is either free or have an electron with 
opposite spin direction. This make the hops longer than in the case when all the sites are 
allowed for the hop. The hnite spin-flip time allows the hop to the site with the same spin 
projection as the hopping electron with some probability p{H). The probability p{H) de¬ 
pends on magnetic held and changes the concentration of sites included in the percolation. 
The concentration of sites is connected to the critical hopping length. The conductivity 
depends exponentially on this distance. Therefore the dependence of the probability p{H) 
on magnetic held leads to the exponentially-strong magnetoresistance. 
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The model used in Ref. 


■B 


10| does not include Coulomb interaction and any of the higher 


order terms (intersite exchange interaction or under-barrier scattering). Also the importance 
of non-equilibrium correlated filling numbers was not mentioned in jsj-Q . Therefore the 
kinetic equations Eqs. (120112^ should be applicable to this case. We have shown that the 
correct procedure is to consider averaged on-site density matrix and only then calculate the 
parameters of percolation like the characteristic hopping length. As a result the magnetic 
held can affect the effective resistances only via magnetization Mi (at least for the case 
of small dc current when resistor network approximation is applicable). When Zeeman 
energy can be neglected in comparison with temperature, the on-site magnetization is absent 
Mj = 0. In this case the spin relaxation time does not contribute to the expression for the 
resistances and does not have any effect on the resulting conductivity. 

The magnetoresistance appear only when Zeeman energy becomes comparable with the 
temperature. In this case the stationary state have hnite magnetization on each site aligned 
along the magnetic held. This on-site magnetization inhuences the resistor network in 
accordance with equation flTIl) and leads to the magnetoresistance. However it is not a novel 
ehect (at least for lar ge H ubbard energy U). It is well known from semiconductor physics and 


was hrst reported in 


36| . It is positive and a linear function of the magnetic held for kBT< 


Ez < ^c^bT, where a is the critical exponent of hopping conductivity and Ez = ^j^bqH is 


the Zeeman energy 27|. At higher helds Ez > ick bT this magnetoresistance saturates. In 
the limit of the small magnetic helds Ez < kBT it becomes quadratic oc {Ez/kBTY [87 1. 

We want to note however that our treatment is directly applicable only to the systems 
that are close to equilibrium. The resistor network explicitly assumes expansion over small 
applied voltages. Also the Hartree decoupling applied to get the kinetic equation can be 


strictly justihed only near equilibrium. There is a numerical Monte-Carlo simulation 


12 | 


showing that pjEf) can in principle inhuence the dc conductivity. It is important that the 
simulation 12(] deals with strongly non-equilibrium systems with the voltage applied to a 
single resistor at least three times larger than temperature. We believe that in order to 
describe the results of with analytical theory one should directly include the parameter 
of the non-linearity into considerations, because spin-blocking magnetoresistance discussed 


12l | does not appear for small applied voltages. 


We want to compare this result with the results of recent publication 


considering 


magnetoresistance due to the spin precession in the hyperfine fields for another system (the 
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model in 


does not allow the double-occupation but include the interference of different 


tunneling pathes). It is shown that the effect of spin precession on the d.c. conductivity is 
related to the correlations of site hlling numbers that appear in the non-Ohmic regime. It 
is out of the Hartree approximation and does not appear for small voltages. 

Note that we do not provide the explanation of organic magnetoresistance within our 
theory. We demonstrate that in the limit of the small electric held where the reduction 
of the kinetic equations to the effective resistor network is justihed the magnetoresistance 
appears only when magnetization is hnite. In order to describe the magnetoresistance 121 
that appears without average magnetization it is necessary to go to the limit of strong 
electric held where the site occupation numbers may be correlated and our kinetic equations 
are not justihed. 


IV. SPIN RELAXATION DUE TO THE RANDOM FIELDS 


In order to demonstrate that the kinetic equations Eqs. fl2Utl23p are a useful tool to under¬ 
stand the physics of hopping conduction we apply them to the problem of spin relaxation 
in disordered semiconductors with the hopping transport. 

We consider the simple case of neighbor hopping. It means that we assume the temper¬ 
ature to be larger than the width of the distribution of site energies e*. In this case the 
energy disorder can be neglected and the positional disorder dehne the distribution of hop¬ 
ping rates. We show that even in this simple case there is a number of novel effects that were 
not discussed previously. This effects are related to the exponentially broad distribution of 
the hopping rates. 

Numerical studies of the hopping conduction usually consider a set of sites with random 
energies on a lattice. Thus the spatial disorder is totaly ignored and only energy disorder 
is considered. However at least one case is known when the spatial and the energy disorder 
lead to different results. It is the case of slow relaxation in systems with the hopping 
conduction and the strong Coulomb interaction 
on the positional disorder in the present paper. 


. It gives us additional reason to focus 


We consider a set of identical sites that are randomly distributed in space with the 
hopping probability exponentially decaying with the distance between sites rjj, Wij oc 
Hq exp(—2rjj/a), where a is the localization radius. At each site there is a random hy- 
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perfine magnetic field "Hj with the characteristic scale ("H). Initially the system is in the 
equilibrium. At t = 0 all sites obtain small magnetization M* aligned along z-axis. As far as 
it is small the relaxation of magnetic moments does not influence occupation probabilities 
/o, /i, f 2 - The system is dehned by two parameters: the conventional localization parameter 
na^ and the ratio of the characteristic precession frequency in the hyperhne magnetic held 
to the hopping rates {71)/Wq. In the case of the variable range hopping conductivity (when 
the energy disorder is essential) one should also keep in mind the relation between the char¬ 
acteristic size of random energy, Hubbard energy and the temperature. These parameters 
however are out of the scope of our work. 

Let us write the general equation for the spin relaxation 

^ Tj.AfO’ - T.,My, (33) 

j 


Ty = U-y/^'' + + lh-/i 


(i) 


+ 


W. 


-u 

ji 


-n 


(i) 


(34) 


Here T^- is the rate of spin transition from site i to site j. Note that in general case Tij ^ Tji 
even in equilibrium. It is related to the fact that different sites have different equilibrium 
probability /i to be single-occupied and thus different equilibrium magnetization. However 
in the considered problem the sites are equivalent and therefore the spin transfer rate is 
directly connected with the charge transfer rate 


Tij — Tji Qc Rij 


(- 1 ) 


A. Spin relaxation in a pair of sites. 

The elementary source of the magnetic moment relaxation is a pair of sites with different 
local hyperhne magnetic helds. These sites are connected by the spin transition rate T. The 
magnetization dynamics of these sites is described by the equations: 

^+??ixMi = T(M2-Mi), (35) 

at 

^ + ??2 X Ms = T(Mi - Ms). 
at 

Let us discuss the relaxation in the two limiting cases: when the precession frequency is 
larger than the tunneling rate "Hi s 3> T and in the opposite limit "Hi s T. 
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In the case of the strong magnetic held the moments precess aronnd the local helds. The 
relaxation of their absolute values is governed by the tunneling rate T. Let us assume that 
at t = 0, I Ml I = I M 2 1 = M. It is easy to show that relation |Mi| = IM 2 I holds during the 
relaxation. Therefore we may write only the equation for M 

^ = -T(l - cos(ZMiM 2 ))M. (36) 


Here ZM 1 M 2 is the angle between magnetic moments Mi and M 2 . In the case of fast 
precession it should be averaged over precession period. Assuming that \'Hi\ 7 ^ \'H 2 \ yields: 


dM 

dt 


= -T 1 


Mill M211 


cos{Z.HiH2) M. 


Ml M2 

Here Mi||(M 2 ||) is the component of magnetic moment Mi (M 2 ) aligned along the local held 

ni{H2). 


The relaxation rate is zero for the case when both magnetic held and magnetic moments 
are aligned along the same axis and Mi = M 2 . In all other cases it is hnite and proportional 
to T. 

In the opposite case T 3 > "H in the initial relaxation phase t ~ 1/T the diherence of 
magnetic moments Mi — M 2 relaxes with the fast relaxation rate T. After this phase of 
relaxation the moments are diherent only because of hnite local helds and Mi — M 2 oc "H/T. 
However the average magnetic moment M 4 . = (Mi + M 2)/2 cannot relax with the rate T. 
The equation for M+ up to the terms oc "H^/T has the following form: 


^ + H+xM+-i|W_x['H_xMJ]=0. (38) 

Here = {Hi±H2)/2. 

The magnetic moment M_|_ precess around average magnetic held, as it can be seen from 
the second term of the equation. The third term gives a small alternation to this precession, 
but more importantly, it yields relaxation of the absolute value M+ = |M+|. 

^ ^ / (K--K+)M^+-M+)A 

dt 2 T nlnl n\Mi j + 


The relaxation of the magnetic moment is proportional to 'h?/T. It is the motion suppression 
of relaxation well known from j^, where it is related to the electron dihusion. Here we 
show that for considered problem the dihusion over macroscopic distances is not required 
for this suppression. It appears even when tunneling between two sites is considered. 
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B. System without disorder 




FIG. 1: The relaxation of the magnetic moment of a hopping system without disorder, (a) The 
dependence of logarithm of (|m|) and (m^) on time for different values of random magnetic field, 
(b) The dependence of relaxation rate on average on-site magnetic field. 


A system with a large number of sites with low disorder can be characterized by some 
average spin transition rate T. The spin relaxation in such a system is analogous to the 
spin relaxation in a pair of sites. For the case of fast hopping T S> "H the magnetization 
is aligned along the initial magnetization axis and slowly process around average magnetic 
field (that tends to zero for a macroscopic system). As a results the magnetization relaxes 
with the rate oc 

For a large magnetic field slow tunneling cannot keep site moments out of the precession 
around their local fields. Due to this precession the macroscopic moment of the system is 
decreased by a factor of 3, according to Kubo-Toyabe formula 4l|]. On the other hand the 
local magnetic moments on sites remains. The following relaxation of the magnetic moment 
goes with the rate T that does not depend on the magnetic field. 

In Fig. ([1]) the results of numerical solution of Eq. (I3^ for the cubic lattice are shown. 
The considered system has the following set of parameters Tq = 1, na}/^ = 0.5 and n = 1. 
The spin transition rate between neighboring sites is T^eib = 0.018. In the model we use 
linearized kinetic equation considering the initial magnetization to be small Mj(0) <C 1. In 
Fig. ([T]) we plot relative magnetizations mj(t) = Mj(t)/|Mj(0)|. Naturally at the beginning 
of the simulation for every site |mj| = 1 and all m* are aligned along z-axis. In the linear 
case nij follow the same equations as Mj. 

We track independently averaged absolute value of the on-site magnetic moment {\m\), 
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and averaged z-component of the magnetic moments (that correspond to the macroscopic 
magnetization of the system). It can be seen that dnring the hrst phase of the relaxation 
(|m|) and (m^) become slightly different dne to random precession in magnetic held. However 
random magnetic held does not rednce the ratio (m^)/(|m|) below 1/3. The hopping makes 
this ratio even larger. Then the relaxation of (m^) follows the relaxation of absolnte valnes 
of the on-site moments. The relaxation has exponential form with some relaxation rate. 

In Fig.(IT]D) we show the dependence of relaxation rate on average magnetic held. The 
dependence is qnadratic for small helds and satnrates when oscillation in random helds 
becomes mnch faster than the spin transition. This result agrees with the two regimes of 
slow and fast hopping described in j42|. 



FIG. 2: The relaxation of the magnetic moment of a hopping system without disorder for large 
random magnetic compared with Kubo-Toyabe formula. The numbers correspond to the parameter 
H/To. 


At very large magnetic helds a minimum appears in the time dependence of (m^) in 
agreement with On Fig. [2] the initial part of relaxation of (m^) is shown. When the 
precession in the magnetic held is much faster than the hopping it is more instructive to plot 
relaxation not versus Tot but versus Ht. For very strong magnetic held (or very slow hops) 
the computations agree with Kubo-Toyabe formula (Fig. [2]). The hnite hopping rate smears 
this dependence. We want to note that the relation of precession frequency to the hopping 
rate corresponding to the appearance of minimum in (m^) is really large. In the discussed 
results neighbor spin transition rate is ~ O.OISTq and the relation H/Tq = 0.5 correspond 
to the precession frequency that is ~ 30 times larger than the spin transition rate. 
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C. Random positions of sites 


The system with the positional disorder is characterized by the parameter na^ where n is 
the site concentration and a is the localization length. For large na^ ~ 1 the hopping rates to 
the neighbors are of the same order for all sites. In this case the disorder becomes relatively 
unimportant. Without polaronic effects this system exhibit the metal-insulator transition 
and its conductivity is of the band type. Due to polaron formation the transport mechanism 
in a system with low disorder can still be of the hopping type but most interesting part of the 
physics that is usually associated with hopping conduction is absent. Therefore we expect 
that the system should behave similarly to the hopping systems without disorder discussed 
in previous section. 

For the small values of na^, when the characteristic distance to the nearest neighbor 
is larger than the localization length ~ 3> a, the hopping rates (even without 

energy disorder) have exponentially broad distribution. Transport in that case is usually 
described in terms of the percolation theory. The conduction is governed by the threshold 
hopping rate Fperfc that allows the percolation over the macroscopic distances. Most of the 
current is carried by the infinite cluster of the sites that are connected with hopping rates 
F ^ Fperfc- The density of the infinite cluster tends to zero for na^ —)■ 0. This small portion 
of sites that form the inhnite cluster plays important role in the theory of variable range 
hopping conductivity. However for the problem of neighbor hopping the infinite cluster 
density becomes small only for very small localization length < 0.1 43|. 

For the systems with low disorder there are two possible relation of parameters ("H) < T 
and ("H) > T. The first relation leads to the motion suppression of the relaxation with 
relaxation rate ~ ('H)^/T. The second relation leads to the relaxation rate ~ T. When the 
distribution of the hopping rates Fjj and the spin transfer rates is exponentially broad, 
the natural situation is that this inequality have different sign for different hops. 

In this case one can assume that in terms of the spin relaxation the sites should be ar¬ 
ranged into clusters. Inside the cluster the spin transfer rates T > ("H) and the magnetic 
moments in the cluster have approximately the same direction. Different clusters are con¬ 
nected by a relatively slow hopping rates T < ("H). Thus the relaxation will be governed by 
the critical transfer rates ~ (?{). The relaxation due to the fast transfer rates T S> {T-i) 
is suppressed by the motion and transfer rates that are smaller than precession rate T -C {T-i) 
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are too slow to influence the relaxation significantly. Therefore it is natural to assume that 
the dependence of spin relaxation on the characteristic scale of random magnetic held should 
be oc ("H), or may be oc ("H)", where a < 2. Note that a statement that spin relaxations 
is governed by the pairs of sites with the hopping rate comparable to the spin precession 
frequency appears in ^ although only for the pairs of sites that are well separated from 
the rest of system. 

In order to study spin relaxation we perform extensive numerical calculations based on 
the general kinetic equation ([33]). We consider numerical samples with random uncorrelated 
positions of sites. Each site is ascribed by a random local magnetic held. The spin transition 
rates between sites are exponentially decaying with distance 


Tij = Toexp(-2rij/a). 

We use cutoh at some large distance Tmax- Tmax is chosen in such way that each site is 
connected in average with 20 other sites. This choice of Vmax allows us to consider relatively 
large systems up to 10^ sites. More details about numerical computations are given in 
appendix. 

We start from the case when all the magnetic moments are aligned along axis and 
consider the relaxation of the relative magnetizations mj(t) = Mj(t)/|Mj(0)| in time. We 
consider three diherent values of parameter 0.2, 0.5 and 1. For = 0.2 {na^ = 

8 ■ 10“^) the system is deeply in the localized regime even without formation of polarons. 
Percolation hopping rate Tperfc is four orders of magnitude less than Tq. However even for 
this small na^ the density of inhnite cluster is rather high ~ 0.85, yet the amount of sites 
that are out of the inhnite cluster is signihcant. For the value = 0.5 [na^ = 0.125) 
the hopping distances are still larger than the localization length and exponentially broad 
distribution of hopping rates persist. However the density of inhnite cluster is > 0.95 and 
there is only small amount of sites that are out of this cluster. Finally for = 1 there 

is no exponential distribution of neighbor hopping rates. In this case hopping transport is 
possible only due to formation of polarons. 

For the small localization length = 0.2 the percolative value of the spin transition 

rate Tperfc ~ 1-8 ■ 10“^To. We start the discussion from the case when the characteristic size 
of local hyperhne magnetic helds is slightly below of Tperfc, ("H) = 10“^To. The results of 
calculations are shown in Fig. |3l We averaged the results over 50 numerical samples with 
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different random positions of sites (each sample contains 10 ^ sites). 



FIG. 3: The relaxation of magnetic moment for = 0.2 and {H) = IO^^Tq. (a) the values 

of average absolute value of site magnetic moment (|m|) and average z-projection of magnetic 
moment (m^) at the initial part of the relaxation, (b) the logarithms of (|m|) and (m^) for all the 
computed time, (c) the distribution of absolute values of magnetic moments at time Tot • 10 “^ = 
10, 20, 30 and 40. 


At the initial phase of the relaxation t ^ there is little change in (|m|) while 

{rriz) decreases due to rotation of the magnetic moments in the hyperhne magnetic held. 
However, the rotation in the random magnetic held alone cannot decrease (m^) more that 
to 1/3. Therefore subsequent relaxation of (m^) accompanies the relaxation of (|m|). 

Then for a large time interval l-it < 20, the relaxation is exponential. However at larger 
time the exponential relaxation slows down. To understand this behavior we computed the 
distribution of absolute values of on-site magnetization at diherent Tot (Fig. [3] (c)). This 
distribution becomes exponentially broad for Tot > 20 and the average magnetic moment 
is determined by a small number of sites with the largest magnetic moment. This behavior 
is quite natural for a system with the exponential distribution of the hopping times. 

In Ref. |45| Shklovskii has introduced a concept of ’’traps” that appear in the hopping 
systems. The typical trap is a pore with relatively large radius that can randomly appear 
in a random site distribution. There is one site inside the pore that act as a trap. The 
capture/relese rate of the trap is related to the hopping rate from the site in the pore to 


the outside sites. In Ref. 


45| this concept was applied to explain 1// noise in the hopping 


systems. We argue that the same traps can be responsible for the spin relaxation at large 
times. Naturally, the spin on these traps cannot relax faster than Ttrap, where Ttrap is 
the highest spin transfer rate from the trap to the surroundings. In section IIVEI we show 
that the spin relaxation due to these traps is non-exponential and derive the analytical 
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approximation for this non-exponential relaxation. 

Note that the initial exponential part of the relaxation appears to be self-averaged even 
for one numerical sample consisting of 10^ sites. However for larger times it becomes more 
difficult to obtain averaged results. For Tot > 50 • 10^ we observed oscillation in (m^) 
that are due to the lack of disorder averaging. We believe that this slow suppression of 
these oscillations with averaging over disorder realizations is related to the trap-dependent 
relaxation. The characteristic size of the trap that governs magnetic moment at time t grows 
with time. However the probability to hnd corresponding trap exponentially decreases with 
the trap size leading to the poor averaging of magnetic moment at large times. 




FIG. 4: The relaxation of the magnetic moment for = 0.2 and different magnetic fields 

{'H)/Tq (a) and the dependence of the relaxation rate (in the exponential phase of relaxation) on 
the magnetic field (b). 


Let us now discuss the dependence of the relaxation on the characteristic value of the 
magnetic field. This dependence is shown in Fig. 01 It is seen from this figure that the 
final phase of the relaxation at large t is essentially independent of the magnetic field. The 
rate of the exponential part of the relaxation depends linearly on the magnetic field, as it 
was predicted above on the basis of simplified arguments. The characteristic time of the 
transition from exponential to non-exponential relaxation decreases with increasing field. 
For large magnetic field the relaxation is non-exponential at any time. 

Consider now the other values of the parameter The results of computation of spin 

relaxation for intermediate localization length = 0.5 are shown in Fig. [5l The qualita¬ 
tive picture of the relaxation is similar to the case of small localization length = 0.2. 
The relaxation is exponential for the short times and then follows universal non-exponential 
curve independent from the value of random field. The rate of the exponential relaxation is 
proportional to the magnetic field for a wide range of fields 0.01 < ("H)/Tq < 0.2 (Fig. |5](b)). 
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FIG. 5: (a) Spin relaxation for = 0.5 for different magnetic fields. Numbers on plots 

correspond to (b) The dependence of the relaxation rate (in the exponential regime) on 

the characteristic scale of the random field. 


At fields ('H)/To > 0.2 the linear dependence starts to saturate. However for these helds 
exponential part of the relaxation is rather small. In this part spin decrease is less than one 
order of magnitude. For ('H)/To > 0.5 it is impossible to separate the exponential part of 
the relaxation. The conductive cluster (with the above formal dehnition) consist of ~ 98% 
of sites making the concept of percolative cluster ill-dehned. Therefore we do not discuss 
the difference between relaxation of sites within and outside of the percolative cluster for 
v}/^a = 0.5. 



FIG. 6: (a) Spin relaxation for = 1 for different magnetic fields. Numbers on plots cor¬ 

respond to {'H)/Tq. (b) The dependence of the relaxation rate on the characteristic scale of the 
random field. 


Finally at £g. [6] we present the results of relaxation computation for = 1. They are 
quite similar to the relaxation for hoppnig system on the lattice. However the relaxation in 
the limit of strong random held is slightly non-exponential. 
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D. Spin relaxation and the infinite percolative cluster 


We have shown that the spin relaxation in the systems with the hopping conductivity 
is extremely non-homogenous. Spin on different sites relaxes with different rates. Note 
that significance of different sites in other situations, for example for the electrical current, 
is also different. The current is carried by the infinite percolative cluster. Therefore it is 
interesting to know whether the relaxation on the sites of the infinite cluster is different 
from the relaxation on the sites out of this cluster. In our study we use the following 
formal definition of the conductive cluster (for details see 2^). First we find the percolative 
threshold, i.e. the critical distance rth that allows the percolation across the numerical 
sample with hops over distances r < Vth. Then we consider sites to be connected if the 
distance between them is less or equal than rth + a/2 (it ensures that hopping rate between 
these sites is higher or comparable with the critical percolative rate). Finally we find the 
infinite cluster connected with distances r < rth + a/2. We consider this cluster to be the 
percolative cluster responsible for conduction. The density of this cluster appears to be 
relatively high (i^ 0.85) for = 0.2. 



FIG. 7: (a) the relaxation of the magnetic moment inside and outside of the infinite cluster, (b) 

the contribution of sites out of the cluster to spin relaxation in the cluster 

We perform the averaging of the magnetic moments independently for the sites of the 
infinite cluster and for the sites outside of the infinite cluster and compare the results with 
the magnetic moment averaged over all the sites (Fig. [7] (a)). One can see that initially 
(7'ot < 2.5 • 10^) the relaxation of the whole numerical sample follows the relaxation of 
the infinite cluster. It is natural because most of sites belong to this cluster. However for 
Tot > 2.5 ■ 10^ these relaxations start to deviate one from another and magnetization of the 
whole sample appears to be larger than magnetization of the infinite cluster. In this time 






domain the relaxation is governed by the spin relaxation on the relatively rare sites with the 
slow relaxation rate that are outside of the inhnite cluster. The relaxation of the average 
moment slows down while the relaxation of the magnetic moment on the inhnite cluster is 
still exponential until Tot ~ 4 • 10^. 

Most interesting is that for larger times Tot > 4 • 10^, the relaxation of sites in the 
inhnite cluster also slows down. We believe that the reason for this slowing down of the 
relaxation is the spin transfer from slow relaxing traps to the inhnite cluster which leads 
to the re-magnetization of the cluster. Actually, the traps can not lose magnetization by 
themselves. They slowly transfer the magnetic moment to the sites of the inhnite cluster 
where it relaxes. To prove that we remove from the system all sites that do not belong to 
the inhnite cluster and recalculate the relaxation. We compare the results with the average 
magnetic moment on the cluster when all the sites in the system are present in Fig. [3 (b). 
For times Tot < 4 • 10^ the curves are the same except for the small diherence in relaxation 
time. However for T^t > 4 • 10^ the relaxation of cluster slows down when all sites are 
included into the computation and remains exponential when we exclude all sites outside 
of the cluster. Therefore we conclude that ~ 15% sites that are not included in the inhnite 
cluster govern the magnetization dynamics of the inhnite cluster at sufficiently large times. 

In the discussed situation the diherence between relaxation of the average moment and of 
the moment of the inhnite cluster become important only when the magnetic moment of the 
system becomes small (~ 10“"^ from the initial magnetization). We believe that it is due to 
the fact that in our situation most sites are included into the percolation cluster. However 
the density of the cluster tends to zero for na^ —)■ 0. Also it is known that the density 
of the inhnite cluster is considered to be small in the theory of variable-range hopping 


26l |. We believe that when the density of the percolative cluster is small the diherence 


between the mean magnetization and the magnetization of the inhnite cluster should be 
more pronounced. 


E. Analytical approximation for non-exponential relaxation. 

We have demonstrated that even at the small magnetic helds the relaxation of the mag¬ 
netic moment have slow non-exponential tails. At large values of the random magnetic 
helds these non-exponential regime covers most part of the relaxation. We have argued 
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that at least for large times this non-exponential relaxation can be described in term of 
traps, the sites that are rather far from their neighbors. Let us discuss the physics of this 
non-exponential relaxation. 

Let us consider the spin on some trap. When the magnetic moment is transferred from 
the trap to its neighbors it quickly relaxes due to fast hops outside of the trap. Therefore 
each trap i has the relaxation rate 7 * that is the rate of electron hopping from this trap. 

Naturally 7 ^ is proportional to the exp(—2r^gjj/a), where is the distance between 
the trap and the nearest site. It is important that S> a. The characteristic number of 


sites that are effective neighbors of the trap is ~ 47 ™(r^gj^)^a. For large it becomes 


,(*) 


larger than unity. In this case 7 * can be expressed as 


7. = Tc 
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FIG. 8 ; Spin relaxation for = 0 . 2 , {H) /Tq = 2-10 ^ compared with analytical approximation 
(|i0|) and (jiT]) 


To estimate the relaxation one should integrate the exponent e with the distribution 
function of distances to the nearest neighbor. 

, Tina^ 4r 4r^ \ 47 rnr^ 


(|m(f)|) = / 47 r? 7 ,r^exp (—Tote 


2 + - + ^ 


a 


dr 


(40) 


The integral fjTni) can be approximated in the limit S> a. 


(|m|)(f) = mo exp 


Tina 


In'(Tot) 


( 41 ) 
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Fig. [8]shows the comparison of the numerical simulations of the relaxation of the magnetic 
moment for = 0.2, {1-i) /Tq = 2-10“^ with the approximate formulae fl^0|) and fl^Tj) . The 
comparison with Eq. fl4Up contain no free parameters, comparison with Eg. 041 p contain one 
free parameter rriQ. We conclude that there is at least semi-quantitative agreement between 
the simulation in the non-exponential phase of relaxation and approximate formulae. 

Eq. finp is in agreement with the expression proposed in Ref. 1^ for the case of the spin- 
orbit mechanism of relaxation in the limit of strong spin-orbit coupling. Note that the similar 
express^. fl has more free parameters. Although our theory aud the theory preseuted 
in Ref. [14j| deal with different relaxation mechanisms, in the limit of strong magnetic held 
or the limit of strong spin orbit interaction the relaxation in both cases is governed by the 
distribution of hopping rates. As a result different spin relaxation mechanisms lead to the 
similar non-exponential regimes of the spin relaxation. 


V. CONCLUSION 

We derived the kinetic equations for the hopping transport that take into account the 
electron spin and the possibility of double occupation. In the limit of low voltage they 
are reduced to the generalized Miller-Abrahams resistor network. We have applied the 
kinetic equations to the problem of spin relaxation in the positionaly disordered system with 
neighbor hopping due to interaction with the random on-site hyperhne (or fringe) magnetic 
held. We show that the initial relaxation rate is governed by the critical hops with the 
rates comparable with the rate of the spin precession in a random helds. At large times as 
well as in the case of the large random helds the relaxation becomes non-exponential and is 
related to the relaxation of the spin in the traps. The relaxation is strongly inhomogeneous: 
the relaxation of sites in the conduction cluster dihers substantially from the relaxation of 
the sites outside of this cluster. However in some cases the traps can drastically ahect the 
relaxation of the magnetic moment of the inhnite cluster. 
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VI. APPENDIX 1: KINETIC EQUATION COEFFICIENTS WITH DOUBLE OC¬ 
CUPATION AND SPIN 


Here we present the derivation of matrix elements W^yz in the eqnation 0181) 

^ W,,Si)pfp^P- (42) 

j 

where x, y and z can correspond to one of /o, /i, / 2 , My or M^. The straightforward 
calculation of matrix elements with equation flT^ is rather cumbersome. The two-site density 
matrix is a 16 x 16 matrix with 256 matrix elements. Therefore we introduce a trick to hnd 
elements W^yz- 

Let us note that equation flT2l) should be valid for any density matrices px{i) and Py{i)- 
However there are some special cases when this equation is reduced to a more simple form. 
Consider for example tunneling from a site i that is single-occupied with an electron with 
spin up (at some moment ti) to a site j that is free at this moment. It means that at t = ti, 

= M® = i, /y = i. 


Other elements of p{i) and p{j) are equal to zero. The transition probability in this case is 
Wij described in Eq. flT^ . Therefore one can write: 


<i/h 

dt ^ dt 


lUi, 


dt 


-IT,,. 


The derivatives of other elements of are equal to zero. The same situation is described 
with equation fflSj) at the moment ti as: 



^/i./o,/(u) + ^/i,/o,M,(u)) 


dt 




dff' 

dt 


^/o./o,/(u) + ^/o,/o,M.(u)- 


As a result we obtain three equations for the hopping rates Wxyz{ij)- 

Considering tunneling of the electron with other spin projections from a single-occupied 
site to a free site we obtain a set of equations 


^/i,/o,/(u) ± kT/i,/o,Ma(u) = Wij, 

^M„,/o,/i(u) ± fhM„,/o,M/3(u) = ±SayWij, (43) 

Wfojojiiij) i H/pjp^Ma(u) = ~Wij. 
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Equations have the following solution: 

^/i,/o,/i(u) = ^fijoMc = 0; 

^McfoJl = 0 ) Wm^Jo,Mi 3 = ^ij'i 

^/o,/o,/i(*i) = ^/o,/o,M„ = 0. 

This gives us 20 (from 216) elements W^^y^z- 

Similarly we can describe other “simple” cases: tunneling a form double-occupied site to 
a free one, tunneling from a single-occupied site to another single-occupied one and so on. 
Each case gives us a set of equations. As a result using these equations one can End all the 
hopping rates Wxyz- 

VII. APPENDIX 2: NUMERICAL TECHNICS 

Our computations are based on the rate equation fl3^ considering Ma'^ {t) — the on-site 
magnetic moments projections averaged with time-dependent density matrix. Here index a 
corresponds to a cartesian coordinate and k numerates the sites. 

We start with the equilibrium state (i.e. the equilibrium density matrix) of the system. 
At the time t = 0 small non-equilibrium part is added to the density matrix so that each 
site has higher probability to be single-occupied with electron with spin up than to be single 
occupied with an electron with spin down. It correspond to the appearance of the small 
2 ;-component of the magnetic moment on each site Mi°^(0) = <C 1. The consequent 
relaxation follows the kinetic equation that relate different magnetizations Ma'^ between 
themselves. 

In our simulation we consider the numerical samples containing 10^ sites with random 
(uncorrelated) positions inside a cube with the side L. We apply periodic boundary con¬ 
ditions. The equation fl33|) is considered in a matrix form drui/dt = '^jTijmj. Here 
correspond to the relative magnetizations mj(t) = Mi(t)/Mi{0). Indexes i and j stand for 
both the number of site and cartesian coordinate corresponding to the spin projection {x,y 
or z). The equation is solved with simple iterations mi{t + 6t) = mi{t) + 6tTijmj{t). The step 
6t was different in different numerical experiment but was selected to be at least 10 times 
smaller than the spin precession time in a random magnetic field and the minimal possible 
spin transition time Tq The results were averaged over at least 100 disorder realizations 
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(i.e. different realizations of random fields and in the case of random site distribution — 
different positions of sites). 

To facilitate the calculation we used a cutoff for the hopping length. The hopping rates 
exponentially depend on the hopping length therefore very long hops do not contribute 
signihcantly to the system dynamics. So we consider only the hops which are not longer 
then some critical distance. This distance is selected in such a way that each site have in 
average 20 neighbors which are included into the computation. We have checked that the 
results do not depend on the cut-off distance. 
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